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The single diffractive cross section for heavy quarks production is calculated in next-no-leading 
order (NLO) for nucleus-nucleus collisions. Such processes are expected to occur at the LHC, where 
the nuclei involved are lead at <Js = 5.5 TeV and calcium at y/s = 6.3 TeV. We start using the 
hard diffractive factorization formalism, taking into account a recent experimental parameterization 
for the Pomeron Structure Function (DPDF). Absorptive corrections are accounted by the multiple 
Pomeron contributions considering a gap survival probability, where its theoretical uncertainty for 
nuclear collisions is discussed. We estimate the diffractive ratios for single diffraction process in 
nuclear coherent/incoherent collisions at the LHC. 

PACS numbers: 13.60.Hb, 12.38.Bx, 12.40.Nn, 13.85.Ni, 14.40.Gx 
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I. INTRODUCTION 

Over the last years, a number of high energy hard 
diffractive processes have used to shed light on our knowl- 
edge about the QCD Pomeron. One of them is the 
diffractive heavy quark production. Central diffractive 
heavy quarks are quite important signals of possible 
new physics, for instance in diffractive Higgs production, 
H — » bb and WH — > W bb. The knowledge of the charac- 
teristic final states with charm or bottom quarks is funda- 
mental to learn how to isolate signals from background 
for very specific production mechanisms, and also, se- 
lects specific classes of relevant higher-order corrections. 
Moreover, it can be useful either to probe the nucleon 
structure and improve the knowledge of radiative correc- 
tions in QCD. 

In the Regge theory, diffractive processes are described 
in terms of the exchange of a (soft) Pomeron with vac- 
uum quantum numbers Nevertheless, the nature of 
the Pomeron and its reaction mechanisms is still not com- 
pletely known. Since the first measurements on diffrac- 
tive jets in hadronic collisions, it is usual consider the 
hard diffractive processes using the diffractive factoriza- 
tion formalism where the hard scattering resolves the 
quark and gluon content in the Pomeron [2| (the so called 
QCD Pomeron). Systematic observations of diffractive 
deep inelastic scattering (DDIS) at HERA have increased 
the knowledge about the QCD Pomeron, providing us 
with the diffractive distributions of singlet quarks and 
gluons in the Pomeron as well as the diffractive struc- 
ture function Q. In hadronic collisions, a single diffrac- 
tive event is characterized if one of the colliding hadrons 
emits a Pomeron that scatters off the other hadron. Hard 
diffractive events with a large momentum transfer are 
also characterized by the absence of hadronic energy in 
certain angular regions of the final state phase space (ra- 
pidity gaps). 

Our goal in this work is to estimate the single diffrac- 



tive cross section for heavy quarks production in heavy 
ion collisions at the LHC. Specifically, we have that the 
coherent diffractive production of heavy quarks in AA 
collisions is the process A + A — > A + [LRG] + QQ + X, 
where [LRG] stands for Large Rapidity Gap. This kind 
of process exhibits a stronger dependence on energy and 
atomic number, with the diffractive amplitude propor- 
tional to the square of the inelastic one. Therefore, it can 
serve as a sensitive probe of the low- a; dynamics of the nu- 
clear matter. The incoherent diffractive scattering in AA 
collisions is the process A + A -> A* + [LRG] + QQ + X, 
where A* denotes the excited nucleus that subsequently 
decays into a system of colorless protons, neutrons and 
nuclei debris. It measures fluctuations of the nuclear 
color field. Such processes can be investigated exper- 
imentally since heavy flavours are copiously produced 
at LHC energies. For instance, ALICE experiment [3| 
is very well suited to perform these analyses because 
it exploits both electron, muon and hadronic channels, 
having a large rapidity coverage and access to \ow-px 
region. Moreover, it presents excellent tracking and 
vertexing capabilities as well as complementary particle 
identification (TPC,TRD and TOF). Charm production 
can be measured using exclusive hadronic channels (as 
D° —> Kir) or using semi-exclusive leptonic channels (as 
c — > i+X), whereas bottom can use similar channels, e.g. 
B — > e/fi + X and b — > £ + X. In special, the tracking 
capabilities at very low transverse momenta in conjunc- 
tion with very nice particle identification of ALICE stand 
it out as one of the most promissing experiments in the 
physics program of diffractive and electromagnetic reac- 
tions 0. 

In the present calculations, we start by the hard 
diffractive factorization, where the diffractive cross sec- 
tion is the convolution of the diffractive parton distri- 
bution functions and the corresponding diffractive coeffi- 
cient functions, in a similar way as for the inclusive case. 
However, at high energies there are important contribu- 
tions from unitarization effects to the single- Pomeron ex- 
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change cross section. These absorptive (unitarity) cor- 
rections cause the suppression of any large rapidity gap 
rocess being important for the reliability of predictions 
. The multi-Pomeron contributions depend in general 
on the particular hard process. They can be estimated 
by a factor (which is energy and model dependent) called 
survival probability factor Q- At the Tevatron energy, 
yfs = 1.8 TeV, the suppression for single diffractive pro- 
cesses is of order 0.05-0.2 [8T4llj|. whereas at the LHC 
energy for pp mode, v/i = 14 TeV, the suppression ap- 
pears to be 0.06-0.1 |8Hlll|. There is a large theoretical 
uncertainty for the nuclear version of suppression factors. 
We will try to do an educated guess of the suppression 
factor for coherent and incoherent channels. 

This paper is organized as follows. In the next section 
we summarize the main formulas considered to compute 
the diffractive ratios for the hadroproduction of charm 
and bottom. We also present the procedure to estimate 
the nuclear incoherent and coherent cross sections at the 
LHC energies. In the last section, we show the numerical 
results for the inclusive and diffractive cross sections as 
a function of energy and give predictions for the corre- 
sponding diffractive ratios for pp and A A collisions. Dis- 
cussion on the nuclear dependence of cross sections and 
corresponding suppression factors are addressed as well 
as a comparison to photonuclear reactions in AA colli- 
sions is performed. 
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the annihilating partons, s is the partonic center of mass, 
rriQ is the heavy quark mass, p R is the renormalization 
scale. It is calculated as an expansion in powers of a s 
with k — corresponding to the Born cross section at 
order 0(a 2 s ). The first correction, k = 1, corresponds 
to the NLO cross section at 0(a^). The dimensionless 
functions have the following perturbative expansion 
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To calculate the /y in perturbation theory, both renor- 
malisation and factorisation scale of mass singularities 
must be performed. The subtractions required for renor- 
malisation and factorisation are done at mass scale p. For 
instance, since the gluon fusion domains at high energies, 
the fq q (p) function is given by [l2| 



II. HEAVY-QUARK PRODUCTION 

Let us present the main formulas for the inclusive and 
single diffractive differential cross sections for the pro- 
duction of heavy quarks in proton-proton collisions at 
high energies. For the diffractive process, the calculation 
is based on the Ingelman-Schlein (IS) model for diffrac- 
tive hard scattering Q. Accordingly, we will take into 
account absorption effects by multiplying the diffractive 
cross section by a gap survival probability factor. Two 
processes are responsible to heavy quarks production at 
the Leading Order (LO) in perturbation theory: quark 
annihilation (qq — > QQ), where the pair is always in a 
colour-octet sate, and gluon fusion (gg —> QQ), where 
both colour-singlet and octet are allowed. At Next- To- 
Leading Order (NLO), the qg + qg scattering is also in- 
cluded. 

In the inclusive case, the process is described for par- 
tons of two hadrons, h a and hi,, interacting to produce a 
heavy quark pair, h a +hb — > QQ + X, with center of mass 
energy y/s. At LHC energies, the gluon fusion channel 
dominates over the qq annihilation process and qg scat- 
tering. The corresponding signal is the production of a cc 
or bb. The NLO cross section is obtained by convoluting 
the partonic cross section with the parton distribution 
function (PDF), g(x,p F ), m the proton, where pp is the 
factorization scale. At any order, the partonic cross sec- 
tion may be expressed in terms of dimensionless scaling 
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where e =_jq§- The f^\p) and flj'(p) functions can be 
found in [12j. 

The running of the coupling constant a s is determined 
by the renormalisation group, 
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where as = f— and ri\ f is the number of light flavors, 3 
(4) to charm (bottom). 

The total hadronic cross section for the heavy quark 
production is obtained by convoluting the total partonic 
cross section with the parton distribution functions of the 
initial hadrons [l2j 

e a b{s,m 2 Q ) = Yl j dxi J dx 2fi(xi,n F )fj(x2,fiF) 



x Vtj(s,m 2 Q ,p 2 F ,p 2 R ), 



(6) 



with the sum i,j over all massless partons. Here, are 
the hadron momentum fractions carried by the interact- 
ing partons, f?^) are the parton distribution functions, 
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evaluated at the factorization scale and assumed to be 
equal to the renormalization scale in our calculations. We 
checked the size of NLO corrections by numerically calcu- 
lating the total inclusive cross sections for heavy flavors 
pair production using the MRST 2001 LO and MRST 
2001 NLO set of partons [3 ■ It was verified that NLO 
calculations in pp collisions are about a factor 1.5 greater 
than the LO calculation, showing the importance of the 
corrections in NLO. The cross sections were calculated 
with the following mass and scale parameters: p c = 2m c , 
m c = 1.5 GeV, /i^ = mj, = 4.5 GeV, based on the current 
phenomenology for heavy quark hadroproduction fl5| . 

For the hard diffractive processes, we will consider the 
IS picture Q, where the Pomeron structure (quark and 
gluon content) is probed. To the single diffraction case, 
it consists of three steps: first a hard Pomeron is emit- 
ted from one of the protons in a small squared four- 
momentum transfer \t\. That hadron is detected in the fi- 
nal state, and the remaining hadron scatters off the emit- 
ted Pomeron. Partons from the Pomeron interact with 
partons from the other hadron and finally, heavy quarks 
are produced in the final state, from the point-like QQ by 
the soft gluon radiation. The reaction for heavy quarks 
hadroproduction is p + p — >• p + QQ + X. In this ap- 
proach, the single diffractive cross section is assumed to 
factorise into the total Pomeron-hadron cross section and 
the Pomeron flux factor J2|. The single diffractive event 
may then be written as [16| 



da SD {hh^h + QQ + X) 
dx)p d\t a \ 



fw/a(x\p \ \ t a \) 

(j(w + h^QQ + X) 



(7) 



where ij is the Pomeron kinematical variable, defined 



is the center-of-mass en- 



as x^p — s$ /s a b, where 
ergy in the Pomeron-hadron b system and \/s b — \fs is 
the center-of-mass energy in the hadron Q -hadront, sys- 
tem, with t a denoting the momentum transfer in the 
hadron a vertex. 

To obtain the corresponding expression for hard 
diffractive processes, one assumes that one of the 
hadrons, lets say a hadron a, emits a Pomeron whose 
partons interact with partons of the hadron b. So, the 
parton distribution in Eq. ^ is replaced by the convo- 
lution between a distribution of partons in the Pomeron, 
Pfa/wiP,^ 2 ), and the "emission rate" of Pomerons by 
the hadron, fw/h.(xTp,t), called the Pomeron flux factor. 
Its explicit formulation is described in terms of Regge 
theory, and so, the expression for the single diffractive 
cross section for QQ production is written as [16j 
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For the heavy quarks production in nucleus-nucleus 
collisions, two processes can occur: the first one is a co- 
herent process, which is when a nucleus emits a Pomeron 
and partons of that Pomeron interact with partons from 
the another nucleus. In this case, the calculations are 
carried out in similar way as described above, where the 
proton form-factor, 



Ami - 2.8t 
Fi(t) = .\ . | I - 



Am 2 — t 



t 



0.7 GeV 2 



(9) 



is replaced by the nucleus form-factor which is parame- 
terized as [131 



(10) 



where Ra is the radius of the nucleus A (Ra = 1.2 A 1 / 3 
fm). That is, the nucleus- Pomeron coupling has the form 
aw oc Af3 \FA(t)\, where /?o is the quark- Pomeron cou- 
pling. Here, we assume the Donnachie-Landshoff nota- 
tion for the nucleon-Pomeron coupling, Pntp = 3/?o-Fi(*), 
and then for the nucleus-Pomeron interaction it is used 
the additivity of the total nucleon-nucleon cross section 
and replace /?o by Af3$ and the isoscalar magnetic nu- 
cleon form factor F\ (t) by the elastic nuclear form factor 
FAit). This procedure has been used for a long time 
to compute the Pomeron-Pomeron contribution in heavy 
ion collisions [3 [l9[ . Notice that a linear dependence of 
the nucleus-Pomeron coupling on A is strong. However, 
it is a reliable approximation for the hard process con- 
sidered here. On other hand, the incoherent process is 
characterized by an emission of the Pomeron for one of 
the protons in the nucleus. Studies for the A-dependence 
of the incoherent diffractive scattering of symmetric nu- 
clei allows the cross section to be parameterized as [l7| . 
In the impulse approximation, this process can be de- 
scribed as the interaction on a nucleus with a Pomeron 
belonging to a nucleon embedded in the nuclear medium 
of the remaining nucleus. Thus, we have a roughly A 2 
dependence for the incoherent cross section: 



inc 
® A 



A a N 



(11) 



where a = 1 in our case. In the literature, some authors 
use a = 0.7 — 0.8 for diffractive process [2(| and aN is 
the diffraction dissociation nucleon-nucleon cross section. 
This theoretical parameterization for the nuclear depen- 
dence is based on the Regge approach to strong interac- 
tions at high energies. The energies for nucleus-nucleus 
interactions considered here are yfs = 5.5 TeV for PbPb 
beams and y/s = 6.3 TeV for CaCa collisions. 

The procedure to compute the hard nuclear cross sec- 
tion above is still arbitrary. Despite soft diffraction in 
nuclear targets to be reasonably described using Glauber 
theory, the situation is not clear concerning hard diffrac- 
tion in nuclear collisions. As we are considering the 
Ingelman-Schlein Pomeron, there are two possibilities: 
(a) one defines a new Pomeron flux depending on the 
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atomic number and differential cross section normalized 
by the Pomeron-nucleon cross section, a Wy>1 or (b) one 
defines a sort of Pomeron- nucleus cross section, cr ffA , 
keeping the original Pomeron flux unchanged. These is- 
sues were first addressed in Ref. 12 1| . where an analy- 
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o- i ac(v / s = 14TeV) 


<7diff(^= 14TeV) 


-Rdiff 


cc 


7811 [fib] 


178 [fib] 


2.3 % 


bb 


393 [fib] 


7 [fib] 


1.7 % 



sis of diffractive dissociation of nuclei in proton-nucleus 
and meson-nucleus was presented. Here, we have cho- 
sen the option (a) due to its simplicity of implementa- 
tion. We also assume no gluon shadowing in the inter- 
acting nucleus. The introduction of shadowing in the 
diffractive cross section should reduce it by 60-70 % 
as R g (x ~ 10~ 4 ) s» 0.65 (22[, where roughly speaking 
x ~ rriQ I y/s. The effect in the diffractive ratio is weaker 
and probably increase it as the new inclusive cross section 
is smaller than the minimum bias approximation. 

In the estimates for the cross sections in Eq. (|5J|, 
we consider a standard Pomeron flux from Regge phe- 
nomenology and which is constrained from the experi- 
mental analysis of the diffractive structure function Q. 
For the diffractive gluon distribution in the Pomeron, 
gw{xi, (ip), we will consider the diffractive PDFs ob- 
tained by the HI Collaboration at DESY-HERA |[, 
where the Pomeron structure function has been modeled 
in terms of a light flavor singlet distribution £(x), i. e., 
the u, d and s quarks with their respective anti-quarks. 
Also, it has a gluon distribution g(z), with z being the 
longitudinal momentum fraction of the parton entering 
the hard sub-process with respect to the diffractive ex- 
change. In our numerical calculations, it will be the cuts 
for the integration over x r ^ m = 0.05. 

As a final step in our estimates of the single diffractive 
cross section, we will consider the suppression of the hard 
diffractive cross section by multiplc-Pomeron scattering 
effects (absorptive corrections). This is taken into ac- 
count through a gap survival probability < \S\ 2 >, which 
can be described in terms of screening or absorptive cor- 
rections [23[ . There are intense theoretical investigations 
on this subject in the last years. We quote Ref. [11| for 
a discussion and comparison of theoretical estimations 
for the gap survival probabilities. We notice that it is 
the main theoretical uncertainty in the present calcula- 
tion of diffractive ratios. As a baseline value, we follow 
Ref. [lOj], which considers a two-channel eikonal model 
that embodies pion-loop insertions in the Pomeron tra- 
jectory and high mass diffractive dissociation. For LHC 
energy on pp collisions, one has < |5| 2 >= 0.06. In sin- 
gle channel eikonal models, this factor can reach up to 
0.081 - 0.086 as discussed in Ref. [ll|]. 

Concerning the model dependence, the single chan- 
nel eikonal model considers only elastic re-scatterings, 
whereas for the multi channel one takes into account 
also inelastic diffractive intermediate re-scatterings. The 
available experimental observables which can be com- 
pared to the theoretical predictions of the survival prob- 
ability factor are the hard LRG di-jets data obtained 
in the Tevatron and HERA [l(| [H| as well as diffrac- 
tive hadroproduction of heavy bosons (W^ and Z°) in 
the Tevatron [24j. Here, some discussion for the nuclear 



TABLE I: The inclusive and single diffractive (corrected by 
absorption effects) cross sections in pp collisions at the LHC. 
The corresponding diffractive ratios, -Rdiff , are also presented. 



case is in order. Currently, there is a lack of informa- 
tion on the gap survival probability in nucleus-nucleus 
and proton-nucleus collisions. It is known that it will be 
much smaller than for proton-proton collisions because 
of multiple collisions of projectile nucleons. In Ref. [25| . 
the central exclusive diffraction Higgs production in colli- 
sions with nuclei is considered, where gluon-gluon fusion 
even in proton-nucleus collisions leads to a very small 
cross section, while 77 fusion gives the dominant con- 
tribution. The conclusion is that the value of the sur- 
vival probabilities is negligible small for gluon-gluon fu- 
sion. The calculation gives < l^ 2 > gg ^H= 8 x 10~ 4 for 
proton-ion collisions at the LHC energy and even small ( 
< \S\ 2 > gg ^H~ 8.1 x 10~ 7 ) for ion-ion collisions. 



A. Results and discussion 

Let us present the results for the inclusive and diffrac- 
tive heavy quarks cross sections for hadronic and nucleus- 
nucleus collisions. The calculations for the inclusive and 
diffractive cross sections, as well the diffractive ratios to 
heavy quark production in proton-proton collisions are 
showed at Tab. (P). We take the value < l^ 2 > pp = 0.06 
for the absorption corrections in hadronic collisions at the 
LHC. The partons PDF and scales are mentioned in pre- 
vious section. For the diffractive gluon PDF, we take the 
experimental FIT A (the fully integrated cross section is 
insensitive to a different choice, i.e. FIT B). The main 
theoretical uncertainty in the diffractive ratio is the sur- 
vival probability factor, whereas uncertainties associated 
to factorization/renormalization scale, parton PDFs and 
quark mass are minimized taking a ratio. The present 
results are consistent with a recent estimation computed 
in Ref. [l6[, where a value < IS"! 2 >= 0.09 was consid- 
ered. In that work, a LO version of cross sections was 
considered and then it can be verified that even the NLO 
corrections are absorbed in a ratio. 

As referred before, in nuclear collisions there are not 
calculations of < |5| > 2 for single diffraction processes. 
Estimations of central Higgs production in pA and AA 
collisions give values much lower than 10~ 4 [25[. Thus, 
for sake of illustration we calculated the diffractive cross 
sections and the diffractive ratios in nuclear collisions us- 
ing a theoretical error band for GSP: the upper value is 
obtained using survival probability for pp collisions at 14 
TeV and the lower value is obtained using the estima- 
tion from Ref. ||, < |5| >\ A = 10~ 6 . The results are 
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Incoherent 


CaCa (cc) 


PbPb (cc) 


CaCa (66) 


PbPb (66) 




1.94 mb 


1.68 mb 


0.04 mb 


0.03 mb 


abs 


186 - 0.003 mb 


4356 - 0.07 mb 


3.78 - 6.3 x 10 -5 mb 


85 - 0.001 mb 


flxnc[%] 


40 


38 


20 


19 


i?fnc[%] 


2.4-4 x 10" 5 


2.28 - 3.8 x 10" 5 


1.2 - 2 x 10" 5 


1.14 - 1.9 x 10" 5 



TABLE II: The incoherent cross section per nucleon for calcium and lead without absorptive corrections. The cross section 
including absorptive corrections (crf^c) and the diffractive ratios are also presented (see text). 



shown in Tabs. ((IT]) and for incoherent and coherent 
collisions, respectively. The ratios are obtained from 

ginc g.coh 
-Rinc = ; RcoYi = (12) 

where <7™ c is given by Eq. (fTTj) and a a = A 2 <jn- For the 
coherent case, a c A oh being the c^v cross section with the 
modification in text below Eq. ([TU|) . As a cross check, 
we found the following values for the inclusive total cross 
sections per nucleon: 4.8 (0.2) mb for charm (bottom) 
production in CaCa collisions and 4.3 (0.17) mb for PbPb 
collisions at the LHC. 

In Tab. (JTTJ) , the incoherent cross section per nu- 
cleon, dine/ A 2 , is presented for calcium and lead with- 
out absorptive corrections. The cross section including 
absorptive corrections, is shown as a band: the up- 
per value is obtained using the gap survival probability 
< \S\ 2 > pp (VI = 14TeV) = 0.06 and the lower one 
corresponds to < \S\ 2 >aa— 10~ 6 [25) . The incoherent 
cross section including the absorption corrections for the 
charm case in PbPb collisions has magnitude of dozens 
of micro-barns in the lower band, which is a hopeful sign 
from the experimental point of view. For bottom, the 
situation is similar with a lower bound of order 1 /ib. 
The diffractive ratios are also presented: i?j nc stands for 
the ratio without considering absorption factor, whereas 
iijj c s stands for ratios taking absorption (the band is 
similar as for the diffractive cross section). The ratio 
without absorptive corrections is almost identical to the 
diffractive ratios for the proton-proton case. The rea- 
son is the A-dependence for the incoherent cross section 
crgjjf cx A 2 <Tpp. The ratio is very small but is still larger 
than for double diffractive production (double Pomeron 
exchange), as computed in Ref. [Vq . 

The coherent cross section and ratios are presented in 
Tab. (jlllj) . using the same notation as for the incoherent 
case. It is noticed that the coherent cross section is larger 
than the incoherent one, reaching a factor 2 for PbPb col- 
lisions. This enhancement is translated to the diffractive 
ratios as well. The reason is the ^-dependence for the 
coherent cross section a co h cx A a a p ^ , where a = 7/3. 
The coherent cross section including the absorption cor- 
rections for the charm case in PbPb collisions reaches 
200 /ib in the lower band, which is still hopeful, and for 
bottom one gets a lower bound of order 3 /ib. 

Let us now comment on the ^-dependence of coherent 
and incoherent scattering. For the incoherent case, we 



consider the impulse approximation where a Pomeron is 
emitted from one nucleon embedded in one of the incom- 
ing nucleus, This Pomeron then interacts with the re- 
maining nucleus, where we have assumed no gluon shad- 
owing xgA{x2, Hq) — Axg p (x2, Hq)- Such approxima- 
tion gives a^f cx A 2 <rj? . If we consider pA collisions, 
the the Pomeron would interact with a proton and we 
would get <jf nc cx A(Tp^. In the coherent case, the A de- 
pendence comes from the Pomeron flux for as nucleus, 
fw /a( x tp i t) cx A 2 |F J 4(t)| 2 (the Pomeron is emitted by the 
nucleus as a whole). After integration over t, one gets 
Zip/a fa*) cx A 2 /R 2 A , with R\ = 1A4A 2 / 3 fm 2 . Once 
again supposing no gluon shadowing in the remaining 
nucleus, one has a^ot cx A 7 / 3 . Using the same argu- 
ments it is easy to see that in pA collisions we would find 
^coh ^ 4 ^ 3 a pp (t ne ^-dependence comes only from the 
integrated Pomeron flux). Notice that for pA collisions 
our calculation is in agreement with the Born approxi- 
mation for the interaction with a single nucleon in the 
Glauber-Gribov approach. Moreover, we verify a para- 
metric enhancement by a factor A 1 / 3 of the coherent 
diffractive cross section with respect to the incoherent 
one in both pA and A A collisions. This fact explains the 
distinct values verified in Table (III) compared to Ta- 
ble (II) and the reason for a different enhancement for 
CaCa and PbPb. Such a similar enhancement has been 
recently found in Ref. [27], where the cross section for 
incoherent and coherent diffractive gluon production in 
qqA collisions (considered a prototype of pA scattering at 
the LHC) was computed in the Color Glass Condensate 
framework. 

Concerning the main theoretical uncertainty on the 
present calculation for the diffractive ratios, some words 
of caution are in order. Of course, the error bands in 
Tables (II) and (III) are a naive oversimplification. The 
suppression factors in nuclear collusions are well known 
to depend on atomic number, center of mass energy and 
on the s pec ific produced final state [28j . For instance, 
in Ref. [28| a careful calculation of absorption correc- 
tions for proton-nucleus collisions was carried out consid- 
ering Drell-Yan and heavy quark production (charm and 
bottom). It was shown that the suppression takes dis- 
tinct values for coherent and incoherent diffraction, being 
the absorption on the coherent case one order of mag- 
nitude stronger than the incoherent one. For instance, 
using the single-channel Glauber approach (which is in- 
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Coherent 


CaCa (cc) 


PbPb (cc) 


CaCa (66) 


PbPb (66) 


0"coh/A 2 


2.9 mb 


3.7 mb 


0.05 mb 


0.06 mb 


abs 
"ooh 


277 - 5 x 1(T 3 mb 


9686 - 0.16 mb 


5.2 - 8.6 x 10" 5 mb 


156 - 0.003 mb 


-Rcoh [%] 


60 


86 


27 


35 




3.6-6 x 1CT 5 


5.2 - 8.6 x 10" 5 


1.6 - 2.7 x 10~ 5 


2.1 - 3.5 x 10" 5 



TABLE III: The coherent cross section per nucleon for calcium and lead without absorptive corrections. The cross section 
including absorptive corrections (o"coh) an d the diffractive ratios are also presented (see text). 



sensitive to the parton longitudinal momentum fraction) 
the prediction gives a factor 4 x 10~ 2 for charm/bottom 
production in incoherent scattering for proton-gold colli- 
sions at yfs = 300 GeV (RHIC kinematics) and a factor 
5 x 10~ 3 for the corresponding coherent scattering. As 
the dependence of the absorption factors is mild on en- 
ergy, the estimation presented in Ref. [25[ for the sup- 
pression in exclusive central diffraction in pA collisions 
(< IS*! 2 > gg ^H= 8 x 10~ 4 ) is either reasonable. The 
situation is however unclear for the nucleus-nucleus col- 
lisions. Hopefully, we can try to do an educated guess 
for the absorption factor in single diffraction cross sec- 
tion in heavy ion collisions. In order to do so, we will use 
the procedure presented in Ref. 29], where the central 
diffraction and single diffraction cross sections in nucleus- 
nucleus collisions are computed using the so-called crite- 
rion C (we quote Ref. [29| for further details). The sin- 
gle diffraction (coherent) cross section for AB collisions 
is given by: 

„SD in ( —in\ „in ( in ^.SD\ 

where a l AB is the inelastic AB cross section considered 
as a function of the nucleon-nucleon total cross section a 
(tJp™ and dp® are the inelastic and single diffractive cross 
sections in proton-proton case, respectively). The expres- 
sion in the equation above can be explicitly written for 
pA collisions where the dependence on the inelastic cross 
section, crp n A (cr) = 1 - [1 -aT{b)] A ~ 1 - exp[-AoT(b)], 
is known in Glauber model for fixed impact parameter b. 
Therefore, using such an information and taking cr™ very 
small in Eq. (fl"3)) . one obtains cr™ = A c g x cr™, with 
A cf f — AJ <PbT A (b) exp[-Aa™T A (b)]. The application 
of the formalism to AA collisions turns out to be difficult 
due to the absence of an explicit expression for a AA (a). 
In Ref. [29l | the authors considered optical approxima- 
tion in which <J l AA is given by corresponding formula for 
pA with A — > AB and an effective profile function for two 
colliding nuclei, T AB = J d 2 bT A (b)T B {b - b). The final 
expression for single diffractive (coherent) cross section 
in AA collisions is given by 

<T S J A (^- S -A)=Al s xa^(VS), (14) 
A 2 cS = A 2 J d 2 bT AA (b) exp [-A 2 0% T AA (b)] (15) 

Using Woods-Saxon nuclear densities and conside ring 
the inelastic cross section cr™ (y/s = 6 TeV) = 73 mb [29[ 



the values for the effective atomic number of the colliding 
nuclei are A 2 S = 6.21 for calcium and A 2 S = 9.52 for lead 
nucleus. 

Numerically, the cross sections obtained for the LHC 
using Eqs. (fP4TT51) are the following: a^ Ph = 1.17 (0.02) 
mb and Oc?Ca = 0.84(0.017) mb for charm (bottom). 
In order to obtaining them we have for the proton 
case cr™(v/s = 5.5 TeV) = 120(2) fib and <t%d{>/s = 
6.3 TeV) = 136(2.8) fib for charm (bottom), where the 
single-Pomeron calculation have been corrected by ab- 
sorption factor in proton-proton collisions. Thus, we con- 
sidered < \S\ 2 >* MR = 0.073 (0.07) for = 5.5(6.3) 
TeV, which is obtained using a parametric interpolation 
formula for the KMR survival probability factor in the 
form < |S| 2 >= a/[b + ln( V / s7s u )] with a = 0.308, b = 
—4.42 and sq — 1 GeV 2 . This formula interpolates be- 
tween the (single diffraction) survival probabilities 10 % 
at Tevatron and 6% at the LHC. The single diffractive 
cross section can be also estimated for pA collisions us- 
ing similar procedure. It is obtained CpPh = 

0.76 (0.018) 

mb for charm (bottom), where A e s — 4.39 (29[. Hav- 
ing the corrected values for the diffraction dissociation 
cross section, we are now in conditions to estimate the 
overall suppression factor in coherent case. In order to 
do so, one takes the ratio between the corrected cross 
section and the single-Pomeron calculation show in the 
first row of Table (III). This gives a suppression factor 
of S A = 7.3 x 10- 6 for charm and S A = 7.7 x 10~ 6 for 
bottom production in coherent single diffraction in PbPb 
collisions at the LHC. These values are not so far away 
from the estimation of Ref. (25[, where central diffrac- 
tive Higgs production is considered. Following Ref. [28j, 
the suppression factor for incoherent diffraction should 
be one order of magnitude larger than for the coherent 
case. This fact has an important consequence as it com- 
pensates the smaller incoherent cross section. The en- 
hancement of coherent to incoherent is proportional to 
A 1 / 3 as referred before and this gives a factor around ten 
for lead nucleus. 

Finally, we compare the present calculation with the 
(inclusive) photoproduction of heavy quarks in ultrape- 
ripheral collisions (UPC's) in heavy ion collisions [30] ]. 
The reason is that the final state configuration is simi- 
lar to the coherent diffraction (both nuclei remain intact 
and one rapidity gap). For instance, in Refs. [HI HH 
the inclusive charm and bottom production is computed 
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in coherent heavy ion interaction considering distinct 
theoretical approaches for the heavy quark production. 
It was found <Jp^ Ph ~ 633 — 2079 mb for charm and 
°PbPb = 8-9 — 18 mb for bottom at the LHC, where 
the lower bound corresponds to the result considering 
saturation model and the upper limit corresponds to 
/^-factorization approach. In pA collisions one has 
OpPk = 5—17 mb for charm and cr^p^ = 81 — 155 fib for 
bottom (we quote Ref. [33] for further details). There- 
fore, for heavy ions the photoproduction channel should 
dominate over the single diffractive channel. The situa- 
tion is completely the opposite in the proton-proton case. 
Recently, the inclusive photoproduction of heavy quarks 
was computed for pp collisions at the LHC [3J, [35| us- 
ing distinct saturation models, which provides us with 
the estimations a^ c = 7.54 — 3.66 fib for charm and 
a pp C = 0-16 — 0.05 fib for bottom. These values can 
be compared to Table (I) and it is verified that the sin- 
gle diffractive channel dominates over photoproduction 
channel. The enhancement of UPC's compared to sin- 
gle diffractive channel in AA reactions is easily under- 
stood from their distinct A-dependences: photoproduc- 
tion grows as oc A 3 (Z 2 enhancement from the equiv- 
alent photon flux plus an additional enhancement from 
the photonuclear cross section, a 7 A °c Aa lp ), whereas 
single diffractive cross section gets a factor A 2 S ~ A 1 / 3 
as shown in Eqs. (I14j|15[) . 

In summary, we have presented predictions for diffrac- 
tive heavy flavor production in heavy ion collisions at 
the LHC. The cross sections are large enough and can 
be investigated experimentally. In our calculations, the 
Ingelman-Schlein picture for hard diffraction was consid- 
ered further corrected by absorption corrections given by 
gap survival probability factor. For the Pomeron struc- 
ture function, the HI diffractive parton density functions 
were considered. We investigate the theoretical uncer- 
tainty on the multiple interaction corrections for the nu- 
clear case and addressed the coherent and incoherent 
scatterings. We are aware of the limitations and theo- 
retical incompleteness of such a picture. However, we 
think it is reasonable for a first exploratory study. An 



alternative approach would be consider the nuclear ver- 
sion for the diffractive production of heavy quarks within 
the light-cone dipole approach of Ref. [36| . In that work 
novel leading twist mechanisms of diffractive excitation of 
heavy flavors in hadronic collisions are proposed, which 
broke the factorization leading to higher twist diffrac- 
tion. Returning to the present investigation, the main 
results are the estimations for AA and pA collisions: we 
obtained cr|£ Pb = 1.17 (0.02) mb and a^ h = 0.76 (0.018) 
mb for charm (bottom). As a byproduct, we also esti- 
mate the overall suppression factor in coherent diffrac- 
tion, which reaches Sa — 7 x 10 -6 for heavy quark pro- 
duction in PbPb collisions. The corresponding factor for 
incoherent diffraction would be one order of magnitude 
larger, which it would compensate the parametric en- 
hancement by a factor A 1 / 3 of the coherent diffractive 
cross section with respect to the incoherent one in both 
pA and AA collisions. We verified that the single diffrac- 
tive channel dominates over photoproduction channel in 
proton-proton case, whereas it is one or two orders of 
magnitude smaller than photoproduction in heavy-ion 
collisions. The enhancement of ultraperiphetral collisions 
compared to single diffractive channel is driven by their 
different A- dependences. Notice that the present calcu- 
lation we provided the fully integrated cross sections and 
an accurate study using relevant kinematic cuts (rapidity 
gap separation and transverse momentum spectrum) in 
the LHC is in order. 
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